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ABSTRACT 

A strong external dc magnetic field introduces a basic anisotropy 
into incompressible magnetobydrodynamic turbulence. Tbe modifications that this 
is likely to produce in the properties of the turbulence are explored for the 
high Reynolds number case. The conclusion is reached that the turbulent spec- 
trum splits into two parts : an essentially two-dimensional spectrum with both 

the velocity field and magnetic fluctuations perpendicular to the dc magnetic 
field, and a generally weaker and more nearly isotropic spectrm of Alfven waves. 
A minimal characterisation of the spectral density tensors is given. Similari- 
ties to measm-ements from the Culham-Harwell Zeta pinch device and the UCLA 
Macrotor tokamalc are remarked upon, as are certain implications for the Belcher 
and Davis measurements of magnetobydrodynamic turbulence in the solar wind. 
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I. INTRODUCTION 


We now have in hand the beginnings of a theoiy of high-Reynolds- 

number msgnetohydrodynamic turbulence^”^^ that is at a level of description 

that is as systematic and inclusive as the corresponding theory for fluid t\xr- 
12 

■ bulence. A significant limitation of this body of theory is that it assumes 
a high degree of symmetry in the statistics of the turbulent fields : spatial 

homogeneity, rotational isotropy, and frequently, temporal stationarity and 
mirror (reflection) invariance. It seems clear, however, that much magnetohy- 
drodynamic turbulence will not be so highly symmetric, both in asti'ophysical 
and laboratory situations . 

The rotational isotropy assumption in particular is limiting, because 
many of the most interesting cases involve an externally-imposed dc magnetic 
field which selects a particular direction in space. While it may be reasonable 
to assxune rotational isotropy about the direction of the mean magnetic field, 
it probably is not reasonable to assume it about the other two directions . 

12 

Moreover, while in ordinary fluid mechanics, turbulence isotropizes itself 
at the smaller spatial scales , it is likeo.y that anisotropy in magnetohydrody- 
namics will persist over the full range of scales to which raagnetohydrodynamics 
is applicable. 

The present paper is intended to propose a qualitative picture of 
turbulent, homogeneous, incompressible, magnetohydrodynamic fluctuations in 
the presence of a strong dc magnetic field = B^. Here, strong implies that the 
energy density associated with is large competred to the sum of the fluctuating 
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magnetic energy density and the kinetic energy density associated with the 
fluid motions of the magnetofluid. V/e will he particularly concerned with 
the differences between this case and the purely isotropic case in which there 
is no preferred direction. 

Convincing data on raagnetohydrodynamic turbulence are still rare. 
Probably the best measurements to date are those from the Culham-Harwell Zeta 
toroidal pinch device. Strikingly similar to some of the Zeta results 
are some recent measurements reported for the UCLA Macrotor tokaraak. ’ ' 

There is not very much help yet from numerical simulations: most of the pub- 

lished results are two-dimensional. Some recent three-dimensional isotropic 
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results have begun to be generated by the group at the Observatoire de Nice, ’ 

but are not yet available as of the time of this writing. It appears that it 

is the case that existing raagnetohydrodynamic turbulence computations are a?l 

spectral-method (Galerkin approximation) computations, assuming rectangular 

periodic boundary conditions and no externally-imposed magnetic field, and do 

not fully address the issues addressed here. Previous analytical calculations 

have addressed the subject of the anisotropies introduced by external dc mag- 
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netic fields for the case of low magnetic Reynolds numbers. ’ ’ 

We rely on a mixture of perturbation theory, model calculations, and 
physically-informed guesswork to arrive at a picture of incompressible magneto- 
hydrodynamic turbxilence, homogeneous but anisotropic due to the presence of a 
strong dc magnetic field The picture must be regarded as conjectural until 

more experiments are done. It is, however, suggestive and compatible with what 

„ , 13,1^,15 , ,, . i6,1T 

is known from Zeta and Macrotor. 


¥e may start from the observation that if the external is strong 
compared to the mean fluctuating field, more energy is required to bend and 
stretch field lines than to translate them, particularly at the larger spatial 
scales, This leads to the feeling, confirmed by the asymptotic analysis of 
Sec. II, that a strong external suppressec spatial variations of the mag- 
netic and velocity fields along the direction. Also, the magnetic fluc- 
tuations and electric field fluctuations are primarily in a direction perpen- 
dicular to The spectrum is conjectured, in Sec. Ill, to consist of a 
highly anisotropic part of the geometry just described pliis a smaller, more 
nearly isotropic part which can properly be described as waves. The tendency 
of a strong magnetic field to enforce two-dimensionality renders several. 
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recent two-dimensional calculations ’ and compUs. ^.Hions'^’ ’■ more generally 
applicable than they might otherwise be. Experimental comparisons ai’e also 
remarked upon, in Sec. III. 

Tui’bulence which fulfills all four synmietries emnnerated in the 
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first paragraph requires only one scalar function to characterize it. When 
the syramrtries are relaxed, a more elaborate characterization is required. 
Section IV is devoted to presenting a framework in which homogeneous but aniso- 
tropic and non-mirror~symmetric turbulence may be characterised. We restrict 
oui’selves in this section to the case in which no net electric current flows 
through the raagnetofluid. Wet cxirrents preclude^ spatial homogeneity and 
periodic boundary conditions, and are thus an additional complication in many 
interesting cases; we postpone some considerations associated with situations 
with net dc current fluxes to a later paper. Section V summarises our conclu- 


II. THE LIMIT OF LARGE 

u 


We consider an incompressible magnetofluid of uniform density. The 
magnetic field is + B(x,t) and the velocity field is y(x,t). Ttie electric 
current density is ^ x B(x,t). Both y and 3 are solenoidal vectors. 

We write the dynamical equations in a common set of dimensionless variables: 

3B 

~ = B » Vv - V • VB + 3« • Vv (l) 

at " ^ 

9v 

s „ . Vv + J X 3 + j X 3(3 - 7p„ 


s _ Y . Vv + 3 • 73 - Vu t *73. 


( 2 ) 


Dissipative terms involviug viscosity and resistivity have been omitted for 

convenience, but can easily be reinserted, p is the mechanical part of the 
’ m 
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nressvu'e. o = n(x,t) is the total pressure, p ulus (5-. + 3) /2. p is de- 
termined by taking the divergence of Eq. (2) and using 7 • y = 0 to get the 
Poisson equation 


7^p = 7 


v3 - y • 7y] . 




p is the solution to Eq. (3) subject to whatever boui'.dary conditions apply. 
Rote that for rectangxolar periodic geometry, p is a ojaadratic functional of 
3 and v. The ambient .magnetic field is asstimed to be spatially uniform, te.m- 
•corally constant, and in the s-direction: 3,^ = 3, e„. If B,- has a spatial 
variation, it is assumed slow compared to the ether length scales of interes-t 
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2 2 

If we assume that the mean plasma energy density, <v > + <3 >, is 
2 

small compared to Bq , we may try approximating Eq.s. (1) and (2) by their 
linearized versions: 


3B 

mt 

3t 



Vv 


( 4 ) 


3v 

^ S • "73 . (5) 

at 


The general solution t*i Sets. (3) and (5) is a linear superposition 
of Alfven waves which have both B and y varying as e:cp i(k • x + , where 
u), = k • The eauations (l) and (2) become intrinsically nonlinear, how- 
ever whenever k becomes nearly perpendicular to The neglected nonlinear 
terms become lar.^'c-r than the linear ones . As k becomes more nearly perpen- 
dicular to the time scales associated with the evolution are no longer 




- (kB^)"**, but are determined by the nonlinear processes associated with the 
c.uadratic terms in Eq.s, (l) and (2). llormal modes can tell us nothing about 
these nearly perpendicvi3.ar motions , whose time scales remain finite even in 
the limit of infinite Ilote that for strictly perpendicular spatial varia- 

tion, the terms involving Bq drop out of Sq,s . (l) and (2) altogether. 

A somewhat more formal demonstration of the two-dimensionality of 
the dynamics can be given simply by assuming a well-behaved perturbation serie 
solution to 3g.s. (l) and (2) as 3^ -?■«>. res'olt is sufficiently interest 

ing that it: may perhaps be worrh demonstrating in detail. We introduce a 
formal expansion parameter into the 3^ terms of Eqs. (l) and (2) and look 
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at the first few orders in the assumed pertui’hation series v - + 

+ ... and B « + Note that B^°^ is not the 

same as g^. We demand that the time and space derivatives remain of 0(1) as 

Bq gets large. What is of interest is what kind of motions are enforced by 

the very large value of B^. Note that what is being effected is essentially 

an expansion of the equations of motion in powers of 1 /Bq, 

The coefficients of successively higher powers of e are equated to 

zero and v;e require that y, B, and j remain solenoidal at each order in the 

expansion. Bor convenience, we invoke rectangular periodic boundary conditions, 

30 that all dynamical vector fields are representable as Fourier series, 

•'I 

The 0(s ) terms simply give that 


and 


jc - 0 




( 6 ) 


( 7 ) 


Equation (o) implies that the seroth-order current must flow along ; J 
e^. Since V • =0 and 


. - 


= I 4'°^ (k,t) expU k • x), 
^ k “ ~ “ 


( 6 ) 


fche non-vanishing Fourier coefficients must have k vectors perpendicular to . 
Written in component notation, 

= (0, 0, (x,y,t)), only. (9) 


j’rom the Fourier-transformed version of 7 x B 


vO) _ ,(0) 


ana tjie assu-mea 


periodic boundary conditions, we have at once that 


( 10 ) 


» L exp(i K *x) 

JlC A# «# 4«w 

» (B^°^ (x,y,t), B^^^ (x,y,t), 0) 


only, vil)h k ' B^°^(k,t) » 0, 


(0) /,. 


Eauation (T) implies that 3v^^V'az « 0, or that 
^(0) ^ (vj°^(x,y,t), V ^^(x,y,t), v,^°^(x,y,t) ) 

X if 


(11) 


wnere 


3v 




X 0. 


3:c 


(■!0\ 


oy 


The geometry implied h. ?.as. (9), (IQi, ixnd (ll), with ail variahles 
independent o? z and the variable magnetic field perpendicular to nas 
terrad as a eonseauanoe of the assumption of a well-behaved expansion as 


en^ 


S -r io, To get the d'/namical behavior of y^’ ' and we need to go on to 

0 

the next order. 

Eq,uating the coefficients of eP to each other gives 


f '> \ 


33 


( 0 ) 


3t 




3y 
0 “3 


( 1 ) 


-( 0 ) 


3y'‘'" fn\ (q) (l) 

no. " i 'O u •"' t " 

SZ ^ ~ ~ ~ 


X O. - vp 

-O 


(0) 


(13) 


(1^^) 


The first three terms of Sc. (13) are independent of s. ^ntagrauing 

' 1 '* 

Sc. (13) over one periodicity lengsh in s at fixed x, y, t, tne r' ' ter.t .*rops 
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8B 


( 0 ) 

'air 




(15) 


The X and y components of Eq. (15) the induction eijuation for two- 
dimensional magnetohydrodynamlcs.^ They are eq,ui valent to the pointwise con- 

(■0) two-dimensional 


servation of the vector potential a 
motion of a fluid element! 


(■^+ a^^^(x,y,t) = 0 


(IS) 


In Ec. (l6) and hereafter, suhscriptn **1." and "jj" will mean ccmponentjs per- 
pendicular to and parallel to respecti'.'-ely . 

Since Eqs. (13) and (15) hold, we ha^'e 3v^“v'az - 0, or - • 

v^^^(x,y,t) only. The s component of Sq,. (15) gives = 0. There 

are two ways this can he achieved. Either v„^^^ » 0, or v„^ ^ « v^ (a ) 
only, which is a statement that v^^^^ is a constant along a fxeld line ot B 
The X and y components of Eq. (ih) give 


3v^°^ 

''i- 


,,(0). ^ ,/0) ^ j(0) gO) _ 

~x 1- ~ '' - 1 - 


\1} 


\~t ) 


+ B |- 

a 32 -1 X 


Before manipulating Eq. (IT), we can obtain one more useful piece 
of information on p^'^^ hy returning to the s component of Eq. vlh): 
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' r\ \ 


•f Y 


V Y 


V t 


f r\ 


r-. - ' 

l i 


The left hand side of P4. (I8) is independent of z. Integrating it over one 
period in a gives 


3t 


+ V 


(0) 


V * 0 

X 3 


( 19 ) 


vhich says that if a fluid element has a s-velooity to start with, the veld- 

t 

city will res.ain constant; no forces act in the z direction. Because of 

a 

So,. (19) » Esi. (18) implies that Sp^^Vsz » 0, or that * P^^^(jc,y,t) only. 

Ue now use the fact that all the terms in iJ<x* (IT) except those 
involving have Been shown to he s-independent . Applying B/Sz to Sq^. (17) 


gives 


^ s + 


Ss 




( 20 ) 


Since 93^^"^V9z + So. (20) implies 


c„(l) 


7 7 
XX ~X 


^u; ^ _~x 


0. 


(21) 


Cm 


,( 1 ) 


Periodic boundary conditions ax*e then invoked again to assert that 3 is 
a supex’-position of linearly independent plane wave contributions of the form 
g,^ e:<p (i k • x) with k • - 0. Eg., (21) implies that 


k k ■ = 0. 


-X-X ~$X 


-!;X 


(22) 


,(1) 


Prom 3c. ^22) we can see that k. and 3, ' are parallel or else both terms in 

'd. 

Eq. (£2) are sex'o. vre must examine these possibilities separately. If both 

texmis and k == 0, then = 0. If the individual terms of Sa* (2'2) are 

s 


not zero, then 


II 


(k^ + k^) k • “ 0 


(23) 


whieh, since 0, implies - 0, ^:mvavy to hj-pothesis. There- 

fore both terms in Sq, (22) vanish separ&telv. For the first term to vanish, 


either » 0 or « 0. If |$ » 0, * 0 and is purely in the 

•** •** “j, Y 

z-direetion and is a function only of s. This cannot he C|^ is solenoidal), 


so k, • * 0 is the only possibility left. ®‘.is implies k •« *1, so that 

can only he a function of n,y, and t. Heturning to (IT), the only 
terms which survive the above conclusions are 


ov, 

-a.,., . 


* .. v” 


X ^ 


3t 


T C24) 

r 


where S is still formally a preasure. It is not necessary 

“ .( 0 ) 


( 1 ) 


(0) , -(0) 


to know ^ to detar.mine » P 




* j 


The Poisson eopaatisn for p^ 


n (3) 

which results from talcing the divercenec ©f S-i. C2-) and using ■' * 0 

determines p' in the usual way. 

3*iuatloa (2^) and she transverse compenenta of 2q. 

(Q) 


"i- 


_ 3 (-).7 V 

*'1 irJ. 


\ J s f-» 


CD 


3t "l ’Iyj_ ~i. Ju 



K J 

system ’ 

whose turb’ilence 
et oj. , ’ 

Wt % 


by Orssag and Tang," and by ? 0 U‘ 4 uet,' from different vlawpoiats, for the case 


\ w J — fS 




j 


The procedure ciui be Iterated to higher order in Ej but the amuaing 
collapae of the geometry to two dimenaional aagnetohydrodynamica does not 
change* 

It should be noted that, though p is an 0(1) (quantity, the Mechanical 

(Ba+B)^ / , -1, 

pressure Pjjj * P r-p*-*. contains formally both 0(e ) and 0(e ) terms. The 

• £« 

first of these makes no contribution to Vp in Ec. (2), but the second Muet be 

in 

-1 

checked in detail to see that there is no Q(e ) contribution , for consistency's 
sake. This in not difficult to do, 



III. QUALITATIVE PICTUBE OP TIJB GPECTEUM 


For turbulence which oceyi spatial homogeneity {i.e., the statistics 
ol* the fluctuating field are translation^invariant) , the most basic variables 
are the covariance tensors of the fields and their Fourier transforms, the 
spectral densities. For example, for the fluctuating part of the magnetic 
field, B, the covariance 

<BCx,t)B(x+r,t+T)> « /dh S(k,T) exp(il$*p), (26) 

where the spectral density tensor G(k,T) may also depend upon the time t in 

-V 

non-stoady situations, is the most physically revealing fiuaatity to consider. 
The trace of S(k,0) meaatires the amount of magnetic energy per unit wavenumber 

jj *» 

space, and various moments of it give such covariances as and <A*B> where 

4 is the vector potential for which 3 “ V x A. Comparable spectral density 
tensors exist for the correlation matrices <w> and <vB>. 

The physics of turbulence is most oatisfaotortly dloeaased in terms 
of the spectral densities, such as S(k,0). These are always the result of a 
balance among three competitive processes: (l) injection of excitations due 

to whatever la driving the turbulexice *, (2) dissipation due to viscosity, 
resistivity, or other decay processes at high wavenumbers 5 and (3) modal trans- 
fer, due to the nonlinear interactions , between one spatial mode and another. 
The range of injection or driving mechanisms is even greater for magnetohydro- 
dynamic turbulence than for Uavier-Stokes fluids , and the various possibilities 
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for modal transfer are also greater; so, depending upon the situation, more 
than one q,ualitative character for S may he expected. Nevertheless, the results 

I 

of Section II suggest a crude qualitative picture of the spectrum of the fluc- 
tuations of homogeneous magnetohydrodynamic turbulence with an anisotropy in- 
troduced by a large value of In Section IV we give a minimal characteriza- 
tion of this field without, however, being able to provide theoretically- 
derived expressions for its elements. 

Ne suggest that a typical turbulent B spectrum will consist of two 
parts. First, the greater part of the energy in k space will reside near the 

plane k-B^ = 0, will involve polarizations such that = 0, and will ap- 

U 5 6 7 Q 

proximate the conditions of two-dimensional ’ ’ ’ ' magnetohydrodynamic tur- 
bulence. It will be non-oscillatory , with time scales which are determined by 
the degree of nonlinearity’ with which the fields are excited. A second part 
of the spectrum will ue more nearly isotropic and can be properly called Alfven 
waves: their time scales will be, predominantly, their length scales divided 

by the Alfven speed constructed from B^. As B^ gets large, they will be 
separated in frequency from the quasl-tvo-dimensional part of the spectrum 
and wi?il have a slower transfer rate, determined by the amplitude of the spec- 
trum itself. In the quasi-two-dimensional part of the spectrum, an Alfven- 
wave-like motion can exist , in which the small-scale fluctuations run along 

22 

the local mean field lines provided by the larger spatial scale components . 

The dynamical role of these pseudo-Alfven waves remains uncertain, but they 

h 5 6 

seem to be effective at enforcing equipartition between magnetic and 


kinetic energies at the smeill spatial scales. In any case, their frequency 
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scales, in the present sltuc-tion. do not vary proportionately to Bq, and they 

are to be sharply distinguished from the true Alfven waves being alluded to 

nere, which are basically the three-dimensional solutions of Eq.s. (4) and (5). 

For what we are calling Alfven waves, these two linear equations are a good 

approximation; for what we are calling the two-dimensional part of the spectrum, 

the linear terns in Eqs, (k) and (5) are identically zero, and Eqs, (2h) and 

(S5) are required. Real life is also likely to involve a transition region, 

» 

with small but non~zero, where linear and non-linear effects will be of 
~ “0 

comparable magnitude. The direction of flow of excitations in k space across 
this transition region is one of the major unanswered questions remaining. 

The partition of the excitations between the two types of turbulence, 
two-dimenslenal and Alfven wave, is bound to be situation dependent and will 
depend upon the excitation mechanism for the turbulence. In laboratory experi- 
ments on confined plasmas, the candidates for excitation mechanisms are very 
numerous. literally hundreds of plasma instabilities (growing linear perturba- 
tions about quiescent laminar states) have been catalogued; at a more elementary 
level, the large radial gradients that are maintained in such fundamental 
paraiv^eters as temperature and pressure, (and frequently, magnetic field, density, 
and fluid velocity) loom immediately, to anyone familiar with fluid turbulence, 
as potential drivers for turbulent motions. Because of the relatively rapid 
variation of the mean properties with the transverse coordinates , compared to 
typically slower axial variations, one may well imagine a selective excitation 
of the two dimensional part of the spectrum. For many turbulence-producing 
agents (impressed changes in the boundary conditions or electrical circuitry 


'm 


1 <? 


supporting the system, for exainple) the time scales will he finite and will 
not speed up as the external magnetic field strength is made larger. One may 
reasonably expect some matching between the time scales of the excited turbU'^* 
lehce and the time scales of the processes which drive it. likewise, it is 
reasonable that the more rapidly vai'ying components of the turbulence will 
usually arise parasitic ally, as a result of nonlinear transfer from modes in- 
volved in the low-freq,uency (i.e. « the reciprocal of an Alfven transit time) 

* 

part of the spectrum. In this case, the Alfven wave component will be regarded 
as derived from the two-dimensional magnetohydrodynamic component. 

One other reason for imagining the Alfven wave component to be 
weaker or less energetic than the two-dimensional magnetohydrodynamic component 
ia that resonance condiviona mitigate against a rapid or effective transfer to 
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the Alfven waves and among the Alfven waves. Speaking loosely for a moment, 
two modes with wave numbers k^^, k^, and fre:^uencieo transfer most ef- 

fectively to a third mode with wavemenber k^ and frequency when m.odal 


matching conditions are met: 


fei + ISa = IS3 


‘^1 ^ ^2 " “3 • 


(27a) 

(2Tb) 


For Navier-Stokes turbulence and two-dimensional magnetohydrodynamic turbu- 
lence, all three of Wg, and always zero, so every ti’iad satisfying 

(27a) is always "resonant". On the other hand, in throe dimensions, some of 
the modes involved may be Alfven waves with w = + k’BQ, and (27b) then consid- 
erably restricts the possibilities for transfer. In particular two two- 
dimensional magnetohydi'odynamic Fourier modes cannot combine resonantly 


to feed an Alfven wave with w„ » k*B^ 5 ^ 3. Two Alfven waves with = au and 
kg'BQ = k^'So resonantly drain a two-dimenaional magnetohydrodynamic aode 
with 0)^ “ 0, hut this is a higher-order process, and therefore a slower one, 
than transfer among three two-dimensional magnetohydrodynamic modea. Of course 
these arguments from the "weak turbulence" perspective are less than rigorous , 
and it is likely that higher-order processes will drain two-dimensional magneto- 
hydrodynamic turbulence into the Alfven wave part of the spectrum, Nevertheless, 

f 

there seems to be qualitative reason to regard the Alfven wave coupling to be 
relatively weak. 

To summarize the above picture, we conjecture several features of 
magnetohydrodynamic turbulent fields existing in a strong do magnetic field 
(l) velocity- field and magnetic-field fluctuations are perpendicular to or 
nearly so; (2) the correlation lengths along B^ are much longer than those 
transverse to since the fluctuating components have little variation 
along Bg', (3) since the electric field E = -y x Bq, the electric field fluc- 
tu ions are also largely perpendicular to and (h) on top of the essentially 
two-dimensional magnetohydrodynamic spectrum is superposed a weaker three- 
dimensional Alfven wave spectrum with frequency scales which sca.Te as Bq. 

One may inquire into the extent to which the above predictions are 

borne out by existing measurements . As far as laboratory measurements go , the 

answer appears to be, rather well. Although the measurements performed were 

in neither case exactly what a theorist would have wished for, two rather dif- 

13 ik 15 

ferent sets of measurements on the Zeta toroidal Z pinch ’ ’ and the UCLA 

Macrotor Tokamak^^ substantiate the above picture in several respects. 

First, the correlation lengths in uhe direction of the mean field (toroidal) 
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direction vere nieaaured to ’’ it le«u?t an order of magnitude greater than the 
cQi^relatioa lengths in the two transverse (radial and poloidal) directions, 
for the magnetic field fluctuations, the electric field fluctuations 
(E 5 -V X Bq, so for the ti’anaverse components, these are essentially velocity- 
field measurements ) , and the electrostatic potential fluctuations . Tlius , a 
high degree of two-dimensionality was indicated in all cases. Second, the 
rms transverse magnetic fluctuations were always larger than the toroidal 
fluctuations : hy more thmi an order of magnitude in the Macrotor measurements , 
and by an uncertain factor in the beta exporimento , partially because apparently 
no distinction was made between parallel fluctuations and toroidal fluctuations 
(the poloidal mean field was lax'ge enough in Seta that it makes a difference) . 
Finally, the single-point frequency measuremento for both the magnetic and 
electric field fluctuations shewed frequency spectra which were rather feature- 
less, were well fit by power laws, and fell off so steeply as to be essentially 
sci'o below either the ion gyrofrequency or the reciprocal of the eorrel.ation 
length divided by the toroidal Alfven speed. This last fact indicates the 
surpriaixig result that not only is the higher-frequoney Alfven wave part of 
the spectinun "weak.", as per c’"* conjectures but that it is in effect neaidy 
absent in those two sitUvations. We have no satisfactoiy explanation for this 
absence. In both sets of measureiaents , the received fluctuating signals were 
filtered of their first several kiloherts before any statiaticval processing 
was carried out, making an assessment of the absolute levels of the rms fluc- 
tuations impossible, since the frequency spectra rise steeply towards nero 
frequency. 


As fai* as space physics measurements are concerned, only some of the 

above conjectures are testable, with data currently available. The best 

measurements available are for solar wind magnetometer data; these are well 
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exemplified by the work of Belcher and Davis and are summarised in the review 
2k 

of Barnes . Some ninety percent of the fluctuating magnetic energy is asao- 

■ P3 

ciated with fluctuations perpendicular to the mean field, but the one 

25 

measiirement of parallel versus perpendicular correlation lengths to date in- 
dicates that the perpendicular correlation length is about a factor of two 
greater than the parallel one, (The ensemble chosen was, admittedly, a rather 
specialiaed one.) The solar wind situation io not as well fit by the above 
analysis as the laboratory one, however, since the mean field is comparable 


to the fluctxiation level . 


IV, MINIMAL GiiARAC'ijBEISATION OF SHE SPECTRA 


If the turbulence were homogeneoua , isotropic , etc . , the spectral 
density tensor in Eii. (Oo) would reduce at t » 0 to 


§(k) s s(k,0) « (I - kk) E«(k)/k‘' 


(28) 


and would be characterised by the single scalar function Eg(k) , the inagnotls 
energy spectrum. It is to the mechanical analogue of this spectral function 
that Kolmogoroff similarity argiuiients are usually applied, leading to the well- 
knowii behavior. In case the oyimaotries are suspended, a more elaborate 

net of dapendenceii upon k is necesssa'cy* 

i 

• For S.,(k), the ijth element of g, homogeneity and reality alone give 

XJ ^ 

if „ 


the conditions S, ,{k) » S..(k) 


U 


‘iV 


'■'41 


(-k). From this it follows that the real 


parts of the siro even under k •*" ~k, the imaglnurj' parts of the 

are odd under k -k, and the diagonal elements are real. 

This restricts the tuuaber of independent functions involved in 3(k) 
considerably, but there are still several. It is clear that the most economical 
choices are desirable in order to repx'Gsent g(k), since any future theory, such 
as a generalisation of the Kolmogoroff 3imilarit5'--\'ariable argvmienta, will 
probably be done on the elements of |. 

After considerable trial and error, we have determined the most 
economical representation of §(k) to be possible in a k-dopendent set of basis 
vectors. A set of configuration-space elements of the tensor have been given 
by Matthaeus and Smith, who elaborate several poi^^to beyond those considered 
here. In particular , if we choose a set of b,asis vectors e^, e,^, whex-e 
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J 


( 29 ) 


% s k :c |./lk X . 


~ ®2 ^ ® 3 * 


then the minimal repz’esentation of S(k) is 

2 


S(k) 


a,i3»l 


cxS ct 3 


(30) 


■The absenee of contributions to the dyadios in Sq. (30) is due to the con- 
dition that the fields be solenoidal: 


k*S(k] 


•k * 0 , 


(31^ 


Even in this k-viependent coordinate s 5 >'stom, it can easily be shora that 
and A,, are real, while A, ^ - A^,, . Thus there are four independent real scalar 
functions involved in this representation of S(k) and there are no additional 
constraints provided by the requirement that the fields be solenoidal. These 
functions are functions of the scalars k, cos 9 (where 9 is the angle between 
k and 3^)* and li (where 1 is the asiauthal angle between tiie pro, Section of (: 
onto a plane rernendicular to 3.,, and some fined direction in than rlane)‘. Tor 

■“v 

turbulence which is isotropic with respaot to rotations about g.^, there is no 
dependence, and the A are functions only of k and cos 9. The limiting 
two-dimensional spectrum -described in the previous two sections is repi-esentel 
by a dominant A.^ which is sharply peaicei near B ~ r/2. Isotropic turbulence 
with no helicity is represented by A,,j = A^. = 0 and with = A,,, being in- 
depen-dent cf CCS d and b* The presence of magnvti-.c helicity ''A’3> where A is 
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the vector potential fluctuation for which B » 7 x A, is signalled hy non- 
zero imaginary parts of A^g. In fact, the magnetic helicity per unit wave 
number space is 




(i/k) - A^g) 


s 2Im A^g/k . 


(32) 


Applying /dk to the expression (32) gives the magnetic helicity <A’B>. Non- 
helical hut anisotropic turbulence is represented by real but unaq.nal values 
of A,, and A^^, with A, ^ » 0. The magnetic modal enersy spectrum can be 
written as A,^ + Agg. 

The coi'reoponding representation of the scalars in x space is much 

more complicated, because the differential equations which result from the 

condition that the fields be solenoidal have no obvious scluticns. Gsnnections 

exist between the various configuration-space-dependent functions which are 
pv oA . 

highly implicit 


iT»23 0f, 25), 


It must be regarded as sr. open question as to u’hether th.ere exists 

12 

any luiivarsal anisotropic spectrum comparable to the Kolmogcroff spectrum 
for the isctropic case. This is as true for anisctrcpic Navier-Stokes turbu- 
lence as it is for magnetohydrodynamics , and it may be a long time before the 
question is answ'ered definitely. It seems likely that if any such universality 
exists, it will be phrased in terms of statements about A., , A^r,» and ImA, 

A cuantity directly expressible in terms of the A , which has been 
measured in the solar wind is the covariance matrix <B{x')|(x}>. In terms of 


A A 




1 


I 


(33) 


<B,(x)B.(x)> * a, '/dk E A e.'E, 
i ~ y~ ^ ’*ct gal « P J 


where are any two unit hasis vectors in a set of k-indepeudent coordi- 

nates. If we assusne that the turbulence is rotationala^ synsnetric about the 
direction of g^, so that the are functions only of k and cos 0, symmetries 

lead to several cancellations of the integrals in Ea* (33). The end result 
is an e:q5ression for the covariance matrix 


<B(x)B(x)> « A 1 -i- S bS 


0 
A+B. 


C3i^) 


where 1 la the unit dyadic and b - 

^ U 


ia a unit x»-ector along the do 
magnetic field, A and B are numbers, integrala of the Tima it follows 

that for turbulence that: ia isotropic with respect to rotations abcut a fixed 
direction, ths principal axis transformation can be used to diatcnalioe 
the covariance matrix with two principal moments eiual and one pi-inaipal aitis 
along 3.,, 


A matrix of the form (At) was not found by Belchc 




0-5 

Davi s'" 


for the covariance matrix <5S> in the solar wind. Sather, the three principal 
moments typically stood in the ratios 5;^:1, indicating pei'hars a two- 
dimensionality to a first approximation, but also a non-sero iepartu’ce from 
isotropy with respect to rotations about any direction. The dir'^ct me-asuremen' 
of such tensors as <3Cx)B(x}> is probably she first single-point msasurer.ont 


wlfcAWW ww U-wwC4»«w.« d.V^ ^ >0*. ..f v««» w ^ ’ 


counters which perait measurements of the velocity field y(jj) permit measure- 
ment of the covariances matrices <v(x)b(x)> end <v(x)v(x)>. The two-point 
matrices <|(x)B(x+r)> are in general necessary to establish the specti’ol 
densities such as |(k), unless some version of Taylor's "frosen flow" hypothesis 
is applicable* 


V* zmmx 


Externally imposed dc magnetic fields introduce a liftsic anisotropy 
into homogeneous I incomprenaihle 5iagnetohydrodj»namic turbulence . In the limit 
of large Bq, the fluctuation spectnim splits into on essentially two dimensional 
part with magnetic and velocity fluctuations nearly normal to and nearly 
independent of the coordinate along plus a more nearly isotropic, weaker, 
highcr-freiiucncy Alfven-wave part. For mvony imaginable turbulence-driving 
mechanisms, the quasi-two dimensional part may be expected to dominate. Such 

13 lU IS 

a spectrum seems to have characterised the 2ota ’ ’ pinch device, well in 


advance any natiofactcry theory, as well as the I'CLA Macrotor tokmaak 


*1 

1 4. i 


flic opcetral density tenoory of tU»‘ covariance matrices can be 
characterlsQd by at most four independent scalar funetiona which , roughly 
opeoitinb', characterise the turbulence by saying how energetic it is, how aniso- 
tropic, and how helical. It is uifxnown as to what extent limiting forms oxiat 
at high Reynolds mutborc for these functions, as they are believed to 
exist for isotropic yavler-Stokes turbulence. 

The most useful futui'c direction for the subject to take might well 
be more thorough-going meacuremeats of the kind that were carried out in Rata. 
Intuitive pictures, involving perhaps no more thaii dimensional considerations, 
need to be developed for relating measured properties of the fluctuations to 
the A^g functions.*^ Finally, the qualitative effects of the different kinds 
cf spectral shapes upon such properties as transport need to bo assessed. It 
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